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Abstract 

A discrete Frenkel-Kontorova model with two alternate spring con- 
stants k\ and hi is studied. It is found that this model has many 
surprising behaviours. The continuum version of this model is dif- 
ferent from the sine-Gordon equation, the continuum version of the 
standard FK model. More interestingly, it has an unpinned commen- 
surate phase which is translationally invariant. A phase transition 
takes plae on a lattice point in the parameter space. We have also 
predicted another interesting phase transition in which an unlocked 
commensurate phase breaking translational invariance will take place 
as the strength of the substrate potential is increased. 
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Many years ago Frenkel and Kontorova [|IJ proposed a discrete version 
of the sine-Gordon equation, which is now known as the Frenkel-Kontorova 
(FK) model. It has been widely studied as a model for many solid-state sys- 
tems 0, ||, f|, ^ such as crystal dislocations, absorbed epitaxial monolayers 
and incommensurate structures, etc. This simple classical one-dimensional 
model consists of a group of particles connected by springs subject to an 
external sinusoidal potential. It has a rich phase diagram [^, 0]. Its com- 
mensurate phase is separated by infinitesimal gaps of incommensurate struc- 
ture. In a commensurate phase, translational invariance is broken and no 
zero-frequency phonon mode exits. However, an incommensurate phase is 
translational invariant and has a zero-frequency phonon mode. There is a 
"transition by breaking of analyticity" || for an incommensurate phase as the 
strength of substrate potential is increased. However, there is no symmetry 
breaking in this transition. 

Many generalized FK models 0, [|], || 10, 11] have been studied for dif- 
ferent purposes. In this paper, we propose to study the following generalised 
FK model, 
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U = I E f(n, h, h)(X n+1 -X n - f(n, a, b)f - — E cos(—X n ). (1) 

4 n Z7T n C 

Here f(n, I, m) = 1(1 — (— l) n ) — m(l + (—1)™). X n is the position of the nth 
particle, a and b are the equilibrium lengths of the two springs with spring 
constants k\ and &2. They satisfy the balance equation k\a = k?\) that is the 
ground state solution of Eq. (|l|) with V = 0. c is the period of the sinusoidal 
substrate potential. The system consists of particles connected alternately 
by two kinds of springs. Obviously, in the absence of a substrate potential 
the ground state is a dimerizated state, as shown in Fig. 1. In many real 
quasi-one-dimensional systems 0] of solid state physics, such dimerization 
often occurs due to lattice distortion as in the case of Peierls instability. 

Differentiating the potential function (1) with respect to X n , we obtain 
the equations for the equilibrium states: 
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hX n+1 - (ki + k 2 )X n + k 2 X n .! - Vsin^fX n ) = 

(2) 

k 2 X n+2 - {h + k 2 )X n+1 + k x X n - Vsm(^X n+1 ) = 
We rewrite them as 

(X n+ i — 2X n + X n _i) + (2g — l)(X n+ i — X n ) — 
(X n -X n _ 1 )-%sinC-?X n )=0 

(3) 

{X n+2 — 2X n+ i + X n ) + (2g — l)(X n+ 2 — X n+ i) — 
(X n+1 - X n ) - ^- Sin (^-X n+1 ) = 

Here g = k\/k 2 . We introduce fields X for a particle with n = even and K 
for a particle with n = odd. In this way, the continuum version of Eq. (^) 
can be written as, 

(4) 

S + (25- X -l)f -f = f«n(?y) 
They are two coupled nonlinear equations. When g = 1, Eq. (f|) becomes 
the sine-Gordon equation [|]. Intuitively, we expect our model to behave like 
the FK model when k\ — k 2 is small. However, we find this is not always the 
case and we will discuss this point in later sections. 

From previous studies 0, H, |J, we know the locked and unlocked phases 
of the FK model are related to the distribution of the particles in the external 
potential wells. For a commensurate phase ||, the particles are symmetri- 
cally placed about the bottom of a well or the top of a well, i.e., there is 
a symmetric point. It is just because of this symmetric point that a com- 
mensurate phase is pinned to the substrate. To reach another degenerate 
ground state of the system, it must move its symmetric point over the top 
of a well to the other side of the well. It has to overcome an energy bar- 
rier; therefore, the system chooses to be locked to the substrate. In contrast 
to a commensurate phase, an incommensurate phase does not have such a 
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symmetric point. Therefore, an incommensurate phase is unlocked and is 
translationally invariant. In this sense, we can regard an incommensurate 
phase as the phase that breaks the symmetry of a commensurate phase. Ob- 
viously, the nonuniformity of the spring constants in our model also has the 
same effect. It makes all the phases less symmetric. Intuitively an incommen- 
surate phase of our model is more asymmetric than that of the FK model. 
So a stronger substrate potential is expected to lock our model. Besides, 
due to the nonuniformity, a commensurate phase possibly does not have a 
symmetric point because it only can exist in some special cases. This means 
our model could have a new commensurate phase which is translationally 
invariant and unpinned. Now we prove it in a more rigorous manner by a 
perturbative study f|. 

For a small amplitude V of a substrate potential, we can view it as a 
perturbation to a uniform ground state of X n — X n _i = a or b. V is treated 
as the perturbation parameter. Then X n can be rewritten as 

X n = E n + 5 n + 6 (5) 

Here E n is the position of the nth particle when V = 0, 9 is an overall phase 
shift factor and 5 n = is equivalent to V = 0. Putting Eq. (||) into Eq. (0), 
we have 

k 2 5 n - (h + k 2 )5 n -i + Mn-2 - Vsin[^{E n -x + + 6)\ = 
k x 5 n+l - (h + k 2 )5 n + k 2 5 n -i - Vstn[^(E n + 5 n + 6)} = (6) 
k 2 5 n +2 - (h + k 2 )5 n+1 + k±5 n - Vsin[^(E n+1 + 5 n+1 + 9)] = 

Expanding 8 n in orders of V, we have 8 n = 5*{V) + S^(V 2 ) + <^(V 3 ) + •••• To 
first order in V, we get 
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k 2 5 l n - (h + k 2 )5l l _ 1 + k x S\_ 2 - Vsin[^(E n ^ + 6)} = 

hSUi ~ (*i + k 2 )5l + k 2 5] l _ 1 - Vsin[^(E n + 9)} = (7) 

k 2 5 l n+2 - (h + k 2 )5i +1 + kt5 l n - Vsin[f(E n+1 + 6)] = 
Eliminating and we have 

<£ +2 - 2# + <£_ 2 = Aszn[— (£„ + 0)] + Bcos[—(E n + 6)} (8) 

where 

A = fcrfe[ fc i cos (v a ) + k 2 cos{^b) + h + k 2 ] 

B = ^-[k lS in(fa) - k 2 sin{^b)\ 

To obtain these results, we have made use of i? n -i = E n — b and i£ n +i = 
E n + a. The solution to Eq. (|8]) is 



= Asin[?f{E n + g)j + (-l)"Bcos[^(£ n + 
2(cos[f (a + 6)]-l) 
Putting this solution into Eq. (@), we obtain 

c2 _ AiSi7l[^(S n + 0)] + (-l^lCOsff (£„ ^ 



(9) 



8(cos [^( a + & )]_l)( cos [^( a + & )]_l) ^ ^ 



with 



A i = rot 5 ( fc i sm (v a ) - k 2 sin(^b)) + A(h + k 2 + hcos^a) + k 2 cos(^b))} 
fi i = iSrMkiSin^a) - k 2 sin(^b)) + + k 2 - k lC os(^a) - k 2 cos(^b))] 

It is not difficult to know to order V n the denominator of 5™ is the multi- 
plication of (cos[^(a + b)} - 1), (cos[^(a + b)} - 1), (cosf^f (a + b)] - 1). 
We immediately know that 5 n is always divergent for any rational (a + 6) / c. 
Unlike the standard FK model, this divergence can simply be eliminated by 
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fixing the phase factor to or c because the numerator of 8™ is a sum of a 
sine and a cosine functions. It seems that perturbation theory || would not 
work; however, this is not the case. In the following, we will show how it 
works perfectly. This divergence only implies there is a new commensurate 
phase which is translationally invariant and unlocked. 

First, we group the solutions 5„_i, S n , 6 n+ \, etc. with two in a group: 
(S n - 2 , S n -i)(i-i),(S n ,S n+ i)i, (8 n+2 ,5 n+3 )(i +1 ), etc. We define a new sequence 
U = 8 n + S n +i- To first order in S n , we have t\ = 5^ + From Eq. @, we 
get 

_ (Acos(f) + Bstn(f))stn[^((a + b)i + 9 + f )] 

cos[f{a + b)]-l 1 ' 

By setting 9 + a/2 = or c, the divergence can be eliminated and tj is locked 
to be identically zero. However, this does not mean the particles are pinned 
to the external potential. It only shows that the two nearest particles move 
in a different direction, i.e., 5\ = —S^ +1 , etc.. Imposing this condition on 
Eq. (0), we get the finite and nonzero solutions 

Vszn[^(E n + 9)] 
K ~ 2(k 1 + k 2 ) {U} 
From Eq. (|T0D, we can extend the above conclusions to second order in 5 n . 
It is obvious that to any order the divergence in S n can be eliminated in this 
way and we can find a finite and nonzero solution like Eq. (|12"D . Summing 
over all these terms, we can write the solution to Eq. (|6]) for any rational 
value of (a + b)/c as 

5 n = Ei{fisin{ 2 -f(E n + 9)} + (-l)»* {Coa [M( Sn + 9 )]}+ 

^m^)SH^{E n + 9)\ W 

Here l(a + b)/c ^ integer, m(a + b)/c = integer, and m(9 + a/2)/c = integer, 
fi is the sum of the solutions to Eq. (§) to all orders that contain the function 
sin[—(E n + 9)]; ti is the sum of solutions to Eq. (||) to all orders that contain 
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the function cos[^(E n + ff)\. For example, to second order in V, fi = 



A/[2(cos[f (o+6)]-l)], h = A 1 /[8(cos[^ (a+b)]-l)(cos[^ (a+6)]-l)] and 
h = B/[2(cos[f (o+6)]-l)], t 2 = B 1 /[8(c OS [^(o+6)]-l)( COS [^(a+6)]-l)]. 



In high orders in V, they are a very complex sum. p m is the coefficient of 
the function of sin[^(E n + 9)] in Eq. (§) just as shown in Eq. (0). It is 
easy to show 



from the conditions m(a + b)/c = integer, and m(9 + o/2)/c = integer. 
So the second sum is independent of the the particle's positions; otherwise, 
this term has different signs for the particles at even and odd sites. This 
shows that all particles with n = even move the same distance in the same 
direction and particles with n = odd move in the opposite direction with the 
same distance. In other words, the two sublattices have a relative parallel 
movement. The distance of the movement is dependent on the overall phase 
factor 9. 9 is determined by the condition m{9 + a/2)/c = integer. It 
is related to the parameter a. However, the parameter a is completely free 
because the equation b = ga does not impose any restriction on it. Therefore, 
the overall phase factor 9 has a lot of freedom in our model. 

Before we go further, let us recall some interesting properties of the stan- 
dard FK model. In the FK model, the perturbative solution for a com- 
mensurate phase can be written as 



Here E n = a n, Iclq/c ^ integer. This solution can be viewed as a special 
case of Eq. (|l"3|) at g — 1. It is an odd function. Therefore, a commensurate 
phase always has a symmetric point. Now let us inspect a more specific case. 
For a /c=M/N, 9/c can take the value n /N with n = 0, N. We have 




(14) 



c 



c 




(15) 



E n + 9 ncio 9 nM + no 



(16) 



c c c iV 
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It is not difficult to find that there is an infinite set of n values that make 
(nM + Uq)/c equal to an integer for any group of M, N and no- This means 
that there is an infinite set of particles that never move with a sinusoidal 
potential V for any commensurate phases and any overall phase factor in the 
standard FK model. They stay at their initial positions for ever no matter 
how 9 is changed. Hence a commensurate phase in the standard FK model 
is always pinned to the substrate and is not translationally invariant. 

Now let us return to our model. The most obvious feature of Eq. (|13|) 
is that there is no obvious symmetry because it is a sum of an odd function 
and an even function. Thus there is a symmetry breaking in our model. A 
commensurate phase in our model does not have any symmetric point. It 
breaks the symmetry of a commensurate phase in the standard FK model. 
This is similar to an incommensurate phase in the FK model. Another inter- 
esting property is that no particles are pinned to their initial positions. This 
is because there is no way to make a sine function and a cosine function to 
be equal to zero at the same value of their variables. Moreover, all particles 
are related through the overall phase factor 9. Its effect cannot be simply 
eliminated as shown in Eq. (|l6l ) or otherwise. Any value of 9 has an effect 
on all the particles. Any change in 9 will make all particles move. There 
is no position that a particle is favoured. We have shown that 9 has much 
freedom. Therefore, a commensurate phase in our model is translationally 
invariant and unpinned. However, the movement revealed in Eq. ( p~3|) is not 
trivial. As we have shown, the second sum in Eq. (|T3|) indicates a relative 
parallel movement of the two sublattices. The first sum also shows the trend 
of two sublattices moving in two different directions because there is a change 
of sign for n = odd and n = even. So the basic movement in our model is 
a relatively sliding between two sublattices. As a result a commensurate 
phase is not expected to have a zero frequency phonon mode for the whole 
system because it requires the whole system to move in the same direction. 
So far we have shown that our model has a new commensurate phase which 
is translationally invariant and unpinned. 
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Except for the new phase, our model also has the usual FK model- 
like commensurate phase for the special rational When k x sin(—a) = 
k2sin(—b),kisin(^-a) = k 2 sin(^-b), etc., the solutions 5^,5%, etc. become 

_ Asin[^(E n + 6)} 
n 2(cos[^(a + b)]-l) 

2 = A lS in[^(E n + 6)} 

n 8 (cos[^(a + b)] - l)(cos[^(a + b)] - 1) 

etc., which are just the same as the solutions f| of the FK model. The only 
way to meet these conditions is a/c and b/c are the two non-zero and non- 
negative solutions of the equation sin(2irX) = 0, i.e., they are integers or half 
integers greater than 1/2. In this case, our model behaves like the standard 
FK model. All the conclusions of the FK model can directly be applied to it. 
These commensurate phases are locked and are not translationally invariant. 

Now we turn to an irrational (a + b)/ c. Like the FK model, the per- 
turbation theory works well and there is no locking phenomenon for any 
irrational (a + b)/c. Therefore, an incommensurate phase of our model 
is translationally invariant and has a zero-frequency phonon mode. For 
a small k, B, B x , ...~0, and A and A x , + cos[^(a + b)}) and 

+ cos\l{a + 6)])(1 + cos[^-{a + &)]), etc.. The solutions (12) and (13) 
become the same as the standard FK model f|. Aside from this, there is no 
way to make our model behave like the FK model for an irrational (a + b)/ c. 

We plot the phase diagram of our model for small V m. a/c versus b/c, as 
shown in Fig. 2. Phases in the parameter space that form a square lattice 
with the lattice constant 1/2 are the locked phases which are not translation- 
ally invariant and do not have a zero- frequency phonon mode. The others in 
the phase space represent unlocked phases which are translationally invari- 
ant. So a phase transition breaking translational symmetry must happen on 
the square lattice in the parameter space (a/c, b/c). 

How this phase diagram changes with the amplitude V of the external 
potential is an interesting question. For an irrational (a + b)/c, no new 
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things are expected to happen: an unlocked incommensurate phase is simply 
locked to a translationally invariant incommensurate phase. However, for an 
unlocked commensurate phase, a new phase transition is expected to occur. 
It will definitely be locked to the external potential with the increase of V. 
It is very possible that an unlocked commensurate phase is locked to a phase 
on the square lattice in the parameter space. If this happens, the phase 
transition would break translational symmetry. 

In summary, we have proposed a new generalised FK model. Its con- 
tinuum version is different from the sine-Gordon equation. Under certain 
conditions, our model is the same as the discrete sine-Gordon equation. Sur- 
prisingly, it occurs at a large difference of (a — b)/c> 1/2 instead of a small 
one for a rational (a + b) / a More interestingly, our model has a translation- 
ally invariant commensurate phase which is unpinned. It has a rich phase 
diagram, which is to be explored in the future. A phase transition break- 
ing translational symmetry take place on a square lattice in the parameter 
space. Furthermore, we predict another new phase transition would happen 
for large V. 
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Figure Captions: 

1. Frenkel-Kontorova model with two spring constants k\ and ki- 

2. Phase diagram a/c vs. b/c. A dot represents a locked phase; the others 
represent unpinned phases. 
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